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ABSTRACT 

The hierarchical growth of mass in the Universe is a pillar of all cold dark matter (CDM) 
models. In this paper we demonstrate that this principle leads to a robust, falsifiable predic- 
tion of the stellar content of groups and clusters, that is testable with current observations and 
is relatively insensitive to the details of baryonic physics or cosmological parameters. Since 
it is difficult to preferentially remove stars from dark-matter dominated systems, when these 
systems merge the fraction of total mass in stars can only increase (via star formation) or re- 
main constant, relative to the fraction in the combined systems prior to the merger. Therefore, 
hierarchical models can put strong constraints on the observed correlation between stellar 
fraction, and total system mass, M^qq. In particular, if this relation is fixed and does not 
evolve with redshift, CDM models predict b = d log /d log A/500 —0.3. This constraint 
can be weakened if the /*-A/5oo relation evolves strongly, but this implies more stars must 
be formed in situ in groups at low redshift. Conservatively requiring that at least half the stars 
in groups were formed by z = 1, the constraint from evolution models is 6 > —0.35. Since 
the most massive clusters (M ^ IO^^Mq) are observed to have f^,^ 0.01, this means that 
groups with AI = 5 x 10^'^ M0 must have /*< 0.03. Recent observations bv lGonzalez et al.l 
Q2007) indicate a much steeper relation, with /*> 0.04 in groups leading to 6 « —0.64. If 
confirmed, this would rule out hierarchical structure formation models: today's clusters could 
not have been built from today's groups, or even from the higher-redshift progenitors of those 
groups. We perform a careful analysis of these and other data to identify the most important 
systematic uncertainties in their measurements. Although correlated uncertainties on stellar 
and total masses might explain the steep observed relation, the data are only consistent with 
theory if the observed group masses are systematically underestimated. 
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1 INTRODUCTION 

An inescapable prediction of all cold dark matter (CDM) models is 
that mass in the Universe, in the form of CDM dominated "haloes'", 
builds up hierarchically, with low-mass systems merging to form 
progressively more massive galaxies a nd clusters of galaxies (e.g. 
iBardeenetalll 19861 : iDavis et alJ[T98^ . Moreover, the rate of this 
mass growth is precis ely determined for a given set of cosmo- 
logical parameters fe.g. lGao"& Whitell2007h . In practice, the com- 
plex and non-linear nat ure of baryonic physic s (particularly c ooling 
and heating processes, [white & Reeslll978l:lwh ite & Freri3 ll99ll : 



heatmg , _ _ . ^ 

ICole et allhoOOl : iBower et alj|2006l : ICroton et alil200 6. and many 
others) means that this is a difficult prediction to test through ob- 
servations of galaxies alone. In fact, observations show that galax- 
ies form in the opposite way, wit h the most massive systems hav- 
ing their stars in place fi rst (e.g. ICowie et alJI 19961 : ljuneau et al.l 
l2005l : |Pozzetti et al.ll2007h . This is not considered a falsification of 



the cold dark matter model, because the effect can be qualitatively 
explained by improving the physical description of baryonic p ro- 
cesses in the models (e.g. lCroton et alj2006l : lBower et al m 
a way that leaves untouched the prediction of hierarchical growth 
in the dark matter component. 

However, an interesting and robust test of the theory can be 
obtained from observations of the stellar fraction of dark-matter 
dominated structures. Unlike gas, it is very difficult to separate stars 
from dark matter, since they are both coUisionless forms of matter 
that interact only via gravity. And while new stars can be formed 
from gas (a process which is very poorly understood) they can only 
be destroyed through normal stellar evolution processes (which are 
quite well understood). The latter effect only removes 10-30 per 
cent of the total stellar mass, with this range reflecting a weak de- 
pendence on star formation history and initial m ass function (e.g. 
Ilungwiert et"al]|200ll : ISruzual & Charlotll2003l) . Therefore, when 
two similar systems merge, the mass fraction in visible stars must 
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be at least as large as the fraction in the combined system prior to 
the merger, and the simplest expectation is that /* will either be 
constant or increase with total system mass. 

Stellar fractions are most reliably measured for galaxy clus- 
ters, where the total mass in dark matter can be determined in 
various, independent ways (e.g. gravitational lensing. X-ray gas, 
or galaxy dynamics). Interestingly, numerous studies have consis- 
tently shown th at ft of clusters and g roups decrease s with increas- 
ing mass (e.g. iHradeckv et al.l [20(M].; Marinoni & HudsonI I2OO2I: 



Eke et alj|2005l: iGirardi et al.ll2002l ; iRamella et al J|2004 iLin et all 
2004lBaloghetalJ 120070 The usual explanation is that clusters 



are built not only from groups, but also from the accretion of low 
mass galaxies, where /, must be very low to e xplain the shal- 
low faint-end slope of the lumino sity function (e.g. IWhite & Fren3 
ll99lUMarinoni & HudsonI l2002h . Another possibility is that low- 
mass groups form a significant number of stars, but only after most 
clusters have been assembled. Either possibility allows theory to 
accommodate a mildly decreasing /, on cluster scales; it is our 
goal in this paper to use conservative constraints on these effects to 
put a robust limit on just how steep this decrease can be. 

An important omission in many of the observational stud- 
ies above, however, has been the contribution from intraclus- 
ter light (ICL), a low-surface brightness distribution of stars in 
groups and clusters that is very difficult to measure. Most stud- 
ies of rich clusters find that the ICL contribution is relatively 
smal l, con tributing less than 3 per cent to the total stellar light 
(e.g.lDurre ll et al. 2002; Feld meier et alj2004lCovone et al.l2003 ; 
I Krick & B ernstein 2 007h. with at most a weak dependence on sys- 
tem mass dZibetti et alj|2005h . Recently, using deep J — band ob- 
servations of 23 nearby svstems. lOonzalez et al.l l2005t) have made 
careful measurements of the ICL component, and come to the sur- 
prising conclusion that both /, and the relative ICL contribution 
depends much more strongly on mass than has been found previ- 
ously, with the ICL actually dominating the total stellar mass in 
groups jOonzalez et al.|[2007l . hereafter GZZ). This has motivated 
us to consider whether or not these observations are able to falsify 
the hierarchical structure growth model. 

We will begin by reanalyz ing the observation al data of GZZ, 
and complementary data from iLin & Moliil ( l2004h in § [H In § [3] 
we use theoretical predictions for the growth of dark matter struc- 
ture to put robust, falsifiable limits on the mass-dependence of /* . 
This prediction is then directly confronted with the observational 
data in §|4l where we also discuss the implications of our findings, 
and the effect of possible biases and uncertainties in the measure- 
ments. Throughout this paper we generally assume a cosmology 
with Q,m ~ 0.3, f^A = 0.7, and Ho ~ 70km s^^. However we 
explicitly consider how our results depend on cosmological param- 
eters, in ij |3.6l 



2 THE STELLAR FRACTION IN LOCAL CLUSTERS 
AND GROUPS 

2.1 Description of the data 

We require an accurate account of the relative stellar content for 
a fair sample of galaxy clusters and groups. We will take most of 
our da ta from two of the best recent surveys, GZZ and Lin & Mohr 
( |2004|) . which are generally complementary in their sources of sys- 
tematic uncertainty. 

GZZ have measured the total light in galaxies and intraclus- 
ter light, for 23 nearby clusters and groups. They use drift-scan 



observation s with careful a t tentio n to flat-fielding, and fit a two- 
component (de^^couleu3 il96lh profile to the brightest cluster 
galaxy (BCG). The outer component, with scale lengths of typi- 
cally a few hundred kpc, is interpreted as the ICL. Stellar masses 
are obtained from the integrated /—band luminosity, assuming 
a mass-to-light ratio of MjLi — 3.6, based on dynamically - 
determined masses for elliptical galaxies l ICapDellari etai]|2006h . 
These masses include a small (about 30%) contribution from the 
dark matter component; thus we adopt a stellar M/L = 2.8 for 
our analysi s, in good agre ement with stellar population models 



assuming a iKroiipal (l200lh initial mass function, as described in 



ICappellari et al.l 1 I2OO6I) . The exact value used does not impact our 
conclusions, which are derived from the trend of the stellar fraction 
with system mass, rather than the normalization. 

For these 23 clusters, the total J— band light from the BCG 
and ICL are well characterized, though the relative contribution 
of each cannot be so robustly determined. The measurement of 
the total galaxy light is made with a statistical subtraction of the 
foreground and background population. The main statistical uncer- 
tainty in these data arises from the total cluster masses, which are 
estimated from the line-of-sight velocity dispersion of the galax- 
ies, a. The relationship between velocity dispersion and dynam- 
ical mass is sensitive to the total potential shape, and velocity 
anisotropy. GZZ partially alleviate this uncertainty by employing 
a calibration between a and X-ray derived m asses, from an inde- 
pendent cluster sample ( IVikhlininet"ai]|200 6^. However, this cali- 
bration sample is very small (13 clusters), so it is not possible to 
determine the scatter about the mean relationship. Moreover, the 
four lowest mass systems in the GZZ sample require extrapolations 
of this calibrating relationship, and may therefore be the most un- 
reliable. 

The second sample we conside r is that of Lin, Mohr & Stan- 
ford (2004) and lLin&Mohj j2004 hereafter LM), who analyze 
if— band observations of 93 X-ray selected clusters, using the 
2 Micron All Sky Survey (2MASS, Jarrett et al. 2000). This has 
the advantages that the stellar mass-to-light ratio {M/L) is only 
weakly dependent on star formation history, and that the total sys- 
tem masses can be determined from the X-ray temperatures. To en- 
sure the stellar masses can be fairly compared with those of GZZ, 
we choose an average M/Lk ~ 0.9; this is consistent with the 
value of A//L/ = 2.8 adopted for the GZZ data if / - A' = 2.0, a 
reasonable number for the early-type galaxies expected to dominate 
these clusters (e.g. IPoggiantill 19971 ; ISmail et alJl200ll ; lKim & Leeil 
I2OO5I ; lEisenhardt et al.ll2007[) . The M/Lk we a dopt is somewhat 
higher than the average value used bv lLin et"al] ( 2004) but, again, 
the absolute value is of little consequence for our analysis. We up- 
date the total mass estimates in LM using accurate A SCA temper- 
ature s available for 63 of their clusters ( Horner 2001 ; Horner et all 
I1999'). These tempera tures are conve rted into MsotQ using the rela- 
tion found by Vikhli nin et al.l l l2006h . where M500 was determined 
from Chandra resolved surface brightness and temperature profiles. 
This typically results in a ~ 20 per cent change to the mass, usu- 
ally in the sense that our new masses are larger. This corresponds 
to a ~ 6 per cent change to 500; therefore we must also correct 
the measurement of total stell ar light within t his rad ius. We simply 
assume the galaxies follow a iNavarro et al. I d 19971) profile, with a 
scale radius = -R500/3, and adjust the total stellar light mea- 
sured by LM accordingly. All of these corrections are small, and 



^ Afsoo is defined as the total mass within a radius 500, such that the aver- 
age density within 500 is 500 times the critical density at the cluster redshift. 
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Figure 1. For the cluster sample of iGonzalez et al.l j2005l) , we show the 
amount of stellar mass in the intracluster Ught, relative to that in the galaxies 
(including the BCG), as a function of total mass Msoo- The solid line shows 
the relation we adopt to correct the data of LM for the intracluster light 
component. 



do not influence our conclusions at afl; however, the more precise 
temperature measurements will result in a more precise mass esti- 
mate which, as we discuss in i; 12.21 is relevant to our interpretation 
of the data. 

Unfortunately, the infrared data of LM are not deep enough to 
measure the ICL contribution directly. To make use of these data 
we need to make an appro ximate correction for t his missing light, 
and we will use the data o f 'G onzalez et"ai] ( |2005[) for this purpose. 
In Figure[T]we show the mass of the ICL, Mid, relative to the mass 
in galaxies (including the BCG), as a function of A/500 • There is a 
correlation (although largely driven by the four poorly calibrated, 
lowest-mass systems) and we find approximately 
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Mbc 



= -0.28 (log A/500 - 12.5) , 



(1) 



shown as the solid line. That is, for the most massive clusters the 
ICL contributes another 20% to the mass observed in galaxies, 
while for the lowest mass systems the total stellar mass is approxi- 
mately doubled by including the ICL. We will apply this correction 
to the data of LM. It should be kept in mind that the ICL and BCG 
light are not robustly separated bvl Gonzalez et al.l ( l2005h . and the 
BCG light in particular is not measured the same way in both stud- 
ies. Nonetheless, the correction is not likely to be grossly incorrect, 
and small differences will not change our main conclusions (see 
further discussion in SHTl. 



2.2 Correlated Uncertainties 

An important consideration in this analysis is proper accounting 
for uncertainties in the measured quantities. In both GZZ and LM, 
the statistical uncertainty on A/*, 500 is small, approximately 10 per 
cent. The dominant statistical uncertainty is in A/500, particularly 
for the GZZ data, where A/500 is derived from velocity dispersions. 
While the number of redshifts per cluster in the GZZ sample is 



generally more than 20, the uncertainties on a are still typically 
~ 10%, which translates to a ~ 30% uncertainty on mass, since 
A/500 oc o"'^- The LM data generally have smaller statistical uncer- 
tainties on the masses, partly due to the fact that the mass depen- 
dence on temperature (A/500CX: T^-^) is weaker than on a. 

However, the error analysis is more complex than this because 
A/*, 500 is the total stellar mass measured within a radius 500, which 
depends on A/500 • Therefore, a statistical overestimate of A/500 will 
also result in an overestimate of A/,, 500, by an amount that depends 
on the radial profile of the stellar mass distribution. For the GZZ 
sample, we use the published r adial profile of the BCG and ICL 
component from lGonzalez et ^ ( 2005). For the galaxy component, 
we only know the total m ass within 500 . We w ill therefore assume 
the galaxy mass follows a iNavarro et alj ( [r997h profile, with a scale 
radius ~ Rsoo/S- For each cluster we can then directly calcu- 
late dA/*,5oo/rfA/5oo and therefore the correlated uncertainty on 
A/., 500, AA/.,5oo = (rfA/*,5oo/dA/5oo) AA/500. For the LM data, 
we do not know the shape of the stellar mass profiles, so we will 
simply adopt the average value of dA/,,5oo/dA/5oo from the GZZ 
data; in any case the error bars for these data are generally smaller 
than the data points in our figures, so this is of no consequence. 

2.3 The stellar fraction in the most massive clusters 

The most massive clusters are the systems for which all matter is 
most reliably accounted for observationally. The hot gas is visible 
as X-ray emission, with a temperature closely related to the grav- 
itational potential, and the total system mass can be directly mea- 
sured via weak-lensing; both of these methods yield masses in good 
agreement with those estimated from galaxy velocity dispersions 
(e.g. lHicks et aL 2006). Moreover, the stellar component is domi- 
nated by old, passive galaxies at the present day, so k-corrections 
and stellar M/L ratios are relatively well determined. 

Considering just the most massive systems in GZZ, those with 
A/500 > 3 X 10" A/0, there are eight clusters with stellar fractions 
ranging from /,= 0.006 to 0.02, and a median of 0.011. In this 
same mass range, LM find that the fraction of mass in galaxies 
spans a range of 0.005 to 0.024, with a median of about 0.0095. 
Including a 20% correction for intracluster light, appropriate for 
these systems (see §[2]l, brings the median stellar fraction to 0.01 1, 
in excellent agreement with the result of GZZ. 

Although these two studies probably provide the most robust 
measurement of the stellar mass fraction in clusters, the results are 
consistent with those of many other studies. For example. E ke et all 
t2004.) find a /3— band mass to light ratio of ~ 350 for the most 
massive systems in the 2dFGRS. Ass uming a typical (but very 
model-dependent) stellar A/ /Lb = 4.5 jFukugita et aljl998h. this 
corre sponds to a stellar fraction in clusters of 0.013. loirardi et al.l 
l2002i) find a higher value of 0.026, assuming the same stellar mass- 
to-light ratio. Most s tudies tend to fin d va lues between these t wo; 
see references within lEke et alj ( |2004|) and lGirardi et alj ( l2002h for 
a good compilation. 

The evidence is therefore good that the stellar fraction in 
clusters is about 1% on average and certainly < 3%. We note 
that this is very similar to the globa l stella r fraction measured 
from large redshift surve ys. lEke et all (l2005l) combined the 2dF- 
GRS dColless et al1l200lh and 2MASS Ijarrett et alJl20 Q0') survej^ 
to fin d an overall stellar fraction of 0.016 (assuming a Kennicutj 
Il983l initial mass function), and t his is consistent with recent re- 
sults from many other studies (e.g. iGallazzi et a"i]|2007l and refer- 
ences within). This small number is already known to put strong 
constraints on the efficiency of galaxy formation when combined 
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Figure 2. The stellar mass fraction /*= M.^soo/A'/soo, including intra- 
cluster light, is shown as a function of total mass M500 ■ The LM data (red 
squares) include a conection for ICL estimated from the GZZ data (circles). 
The two open circles represent clusters A2405 and APMC020, which are 
systems strongly affected by line-of-sight structure. The Icr error bars are 
derived from the published uncertainties on A/500 > and the tilt reflects the 
con'elated uncertainty in A/500 and A/,. 500 /A/500 as described in 12.21 
The horiz ontal, dotted line sh ows the global baryon fraction measured by 
WMAP3 jSoergel et alj2007h . The two solid lines show constant slopes of 
—0.35 and —0.05, for comparison wi th our most conservative theoretical 
lower limit, and the lBower et alj j2006l) model prediction, respectively. 

with measurements of the hot gas mass (e.g. ISalogh et al 
In § [3] we will show that the stellar fraction alone, independent of 
how much hot gas may be present, can be used to test models of 
hierarchical structure formation. 



2.4 Observed Mass-Dependence of /* 

We now consider how /, is observed to depend on system mass. We 
reproduce the data of LM (including our improved mass estimates 
and a correction for the ICL contribution) and GZ^ in Figure [2] 
The error bars represent one standard deviation from the published 
uncertainties on M500, and we include the correlated uncertainty 
on M*,5oo as described in m.2\ However, in this case the direction 
of the correlated errors is dominated by the fact that A/sqo appears 
on both axes. 

Both studies find the stellar fraction decreases with increas- 
ing cluster ma ss, in qualitative agreement with other work (e.g. 
Hradeckv et aU,2() ( )P; Marinoni & Hud son 2002; Ek e et aL 2005; 
Girardi et alji2002l ; iRamella et al.ll2004l ; iBaIogh et alfcoovh . How- 
ever, it is clear that the slope of the relationship is different, and the 
two datasets are therefore inconsistent for the low-mass systems, 
with M < 2 X 1O"A'/0. In fact, GZZ find a remarkably strong 

^ Note that the four lowest mass systems shown here were excluded from 
Fig.l of GZZ, because of uncertainty in the gas mass. However, the gas 
mass is irrelevant for our purposes, so we include these groups here. Two of 
the systems, A2405 and APMC020, have two discrete redshift peaks, and 
A2405 is a clear superposition of two systems. 



trend, with a slope dlog /*/dlog M500 = —0.64, such that the 
least massive systems in their sample have /, > 0.17, in excess o f 
even the total baryon fraction of the Universe ('Spergel et al j2007h . 
Note that the steep trend is not just driven by the last four points, 
however, and the discrepancy exists even if we ignore these groups. 
This is a surprising result and, as we will show in the following sec- 
tion, potentially poses a challenge to hierarchical structure growth 
models. 



3 COLD DARK MATTER PREDICTIONS 
3.1 The progenitor /*-A/50o relation 

The data presented in the previous section suggest that /* decreases 
with increasing mass, above M > 5 x 10^^ Mq. Since CDM the- 
ory predicts a robust connection between structures on these scales, 
we will attempt to secure a prediction for how steep the ft-Msoo 
relationship can be in this regime. 

We start by assuming a logarithmic relationship between /, 
and for systems with M > 5 x 10^^ Af©, at the present 

day. We will then attempt to constrain the slope of this relation, 
b = dlog /,/dlog AfsoQ. The model is normalized so that the 
most massive clusters, with Af5oo= W^^Mq, have /,= 0.01, as 
motivated by observations. To be conservative we will assume that 
haloes with Af < W^^Mq, below the resolution limit of our sim- 
ulations (described in the following subsection), carry no stars, so 
/*= 0. For intermediate masses, lO" < M/Mq < 5 x 10^^, 
where observational constraints are most difficult to acquire, we 
will adopt three ad-hoc models, which span the full range of realis- 
tic behaviour: 

(i) Min:The stellar fraction is assumed to be constant at a low 
value of /*= 0.01. This is the most conservative (non-evolving) 
model, as it will accommodate the steepest slope b while maintain- 
ing consistency with a low /* in massive clusters. 

(ii) Extrap: The ft-Msoo relation is extrapolated to lower 
masses with the same slope b. However we do not allow it to excee d 
the universal baryon fraction of 17.5 per cent te pergel et alj|2007l) . 

(iii) Mirror: /, declines with decreasing mass below 10^^ Mq, 
with slope +|fe|, mirroring the trend at higher masses. Note that 
since this is a logarithmic slope, /* > 0.0 at all masses above our 
resolution limit. 

These three models are shown in Figure|3] with arbitrary slopes (for 
illustration purposes only) b — —0.25, b — —0.40 and b = —0.75 
for the Min, Extrap and Mirror models, respectively. In all cases b 
refers to the slope for M > 5 x 10^^ Af©, and in practice this slope 
is a free parameter that we wish to constrain. 

These models characterise the z — relationship between /» 
and mass. This correlation is likely to evolve with redshift, at a rate 
that depends on the star formation and mass accretion history of 
haloes at a given mass. We therefore consider a generalized model 
of/,: 

log f* = fo (1 + 2)" + fo (1 + zY log Msoo, (2) 

where a and c are free parameters that describe the redshift evolu- 
tion (in normalization and slope, respectively), and /□ is the nor- 
malization of the models at z = 0, fixed so that /,= 0.013 at 
Afsoo = 7 X W^'^Mq. The present-day slope 6 is a free parameter 
for Af > 5 X 10" Af©; below this mass b behaves as described 
above for the Min, Extrap and Mirror models. 
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Figure 3. Examples of our three assumed models for the dependence of /» 
on halo mass at 2 = 0. For M > 5 X 10^^ Mq, the Min model shown here 
(red, solid line) has a slope of —0.25, the Extrap model (blue, long-dashed 
line) has a slope of —0.40, and the Mirror model (green, short-dashed line) 
has a slope of —0.75. The horizontal, dotted line shows the universal haryon 
fraction from WMAP-3 tSpergel et aL2007.) . 

3.2 The simulations 

We obtain merger tr ees generated with the Pinocchio algorithm 
jMonaco etaLll2002h to determine the merger history for an en- 
semble of galaxy clusters. The Pinocchio code has been shown to 
provide results in excellent agreement with full N-body simulations 
(e.g., it reproduces the mass function of simulated dark matter ha- 
los to ~ 10% accuracy). We have verified this agreement, for a 
ACDM universe, by making an explicit comparison of the predic- 
tions of the mean mass growth rate of massive clusters (and the 
scatter about the mean) with that of clusters of sim ilar mass in th e 
publicly available numerical simulation of Spring eTet al.l jlOOSh . 
The Pinocchio code allows us to efficiently change the mass resolu- 
tion (to test the sensitivity of the results to this quantity), to increase 
the simulation box size (in order to generate larger samples of the 
most massive clusters), and to quickly explore other cosmological 
models. For all the runs presented in this paper, we adopt a fixed 
particle mass of 1.35 X 10^°Mq, which is sufficient to follow 
the growth of the massive (> IO^^Mq) systems we are interested 
in. A "resolved" halo corresponds to a system with a minimum of 
10 bound particles. 

We can then compute the final stellar fraction of a halo of 
given mass, based on its merger history, by applying each of the 
models described in §[3T|to its progenitors. As an example, in Fig- 
ure |4] we show the evolution of /, for two massive clusters, using 
a non-evolving {a = c = Q) Mirror model (represented as the 
dashed line) to assign stellar fractions to each merging fragment. 
Each point represents a merger above our resolution limit. At early 
times, /, actually declines, because most of the mass is being ac- 
creted below the resolution limit, where we assume conservatively 
that /,= 0. It is evident from the distance between adjacent points 
that most of the mass is accreted at late times, in low-mass haloes of 
M ~ 1-10 X 10^^ Mq \ hence /* rapidly rises to the corresponding 
value of ~ 0.03, and then remains nearly constant at that level. The 



Figure 4. Two examples of the growth history of massive clusters from our 
ACDM Pinocchio simulations. Each point (small triangles for one cluster, 
and small squares for the other) represents a merger event above our mass 
resolution. The stellar fraction of merging haloes are assumed to follow the 
Mirror model relation shown by the dashed line. For reference, the large 
red, yellow, green and blue points represent the epochs z = 2.0, 1.0, 0.5 
and 0. At early times the cluster is growing primarily from sub-resolution 
accretion; because we conservatively assume this accretion cairies no stars, 
/, slowly decreases. However, most of the mass is accreted at late times, 
through merging with fragments that have M ~ 10^"^ Mq. Therefore, after 
2 ~ 1 /» quickly increases to ~ 0.03, and then remains constant at a value 
higher than the dashed line. The cluster that accretes more of its mass from 
larger progenitors ends up with a lower /« , as expected. 

second cluster trajectory we show is deliberately chosen because 
it is built from more massive fragments and, as expected, it ends 
up with a lower final /* . However, even this cluster has a value of 
/« that is a factor of two larger than the assumed Mirror model. 
Therefore this model is internally inconsistent: as it is not possible 
to preferentially remove stars from the final system, we conclude 
that /, is too high in the low mass systems (i.e., h is too negative) 
to be consistent with these being the progenitors of today's clusters. 
Since we assume /,= below our resolution limit, increasing the 
resolution (i.e., lowering the mass limit) serves to increase /, in the 
final clusters; apart from this, our results are insensitive to changes 
in resolution. 

3.3 Example: CDM predictions for a non-evolving, steep 
/.-A/500 relation (6 = -0.64) 

To demonstrate how our analysis works, we will consider a non- 
evolving model with b — —0.64, which is the observed z = slope 
as measured from the GZZ data. To be conservative we will adopt 
the Min model prescription for haloes with M < 5 x IO'^^A/q; 
this gives the best chance of obtaining low /, in the most massive 
clusters given a realistic growth history predicted by CDM and a 
non-evolving /, -A/500 relationship. In the following discussion, 
the model refers to the assumed /, -A/500 relation that we use to 
specify the stellar content of the progenitors for a given cluster. 
The prediction refers to the f^-M^oo relation that results when a 
given model is applied to the progenitors of an ensemble of clusters 
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Figure 5. Top: An example of the predicted /* for an ensemble of clusters 
in a ACDM Universe, assuming a non-evolving Min model with b = —0.64 
(shown as the red, solid line) for the relationship between /» and M^oq. 
The most massive clusters end up with /« that is much too high, indicating 
an inconsistent model. In contrast, the lower mass clusters lie below the 
solid line, because they are built from systems with /, ^ 0.01. Middle: 
This shows the fraction of stellar mass that has to be added via in situ star 
formation to keep haloes on the assumed Min model (no star formation is 
invoked if /» is greater than this). Bottom: The average redshift at which 
stai's are added to the halo in situ, as a function of final mass. In this case, 
we need to add many stai's at late times in the lowest mass haloes, to keep 
/, as high as 0.05. 



with a CDM-specified merger history. Therefore a self-consistent 
model in the context of CDM is one for which the predicted relation 
agrees with the model relation. 

In the top panel of Figure[5]we show the prediction of the final 
/, resulting from this particular model for an ensemble of clusters. 
The results are somewhat complex, and we see the final clusters 
are grouped into about six "families". Let's consider first the high 
mass end, Af > 5 x 10^*, where the final clusters follow a slope of 
b = —0.64, like the assumed progenitor model, but with discrete 
offsets. The majority of these clusters lie nearly on the model line. 
These are systems which had only one massive 
progenitor, which by construction was assumed to have a stellar 
fraction given by the solid line. A small contribution to the final 
mass comes from lower mass systems (with /*= 0.01) which is 
why the clusters lie slightly below the solid line. The next family of 
points to the right are those clusters that formed from exactly two 
massive progenitors; thus they have similar values of /, compared 
with the first family, but double the final mass. Similarly, the clus- 
ters found farther to the right are made from progressively more 
progenitors with M > 5 x 10^^ Mq. Now considering the whole 
cluster population, the dominant family is actually comprised of 
those clusters with /, < 0.01; these are the systems with zero mas- 
sive progenitors. They have all been built from low mass haloes, 
including a significant contribution from sub-resolution haloes for 
which /,= 0. 

The predicted /, -A/500 correlation is therefore very different 
from the assumed model (solid line). Since we predict a substantial 
number of massive clusters with /« higher than the solid line, this 
model is internally inconsistent; there is no plausible way to remove 
stars from the final system, so we can rule this model out without 
further consideration. 

What about the clusters that are predicted to have a lower 
/, than the assumed model? Unfortunately, these do not represent 
such a clear failure of the model, since one could always assume 
that systems form enough stars in situ, after a merger event, to move 
them up to the appropriate value of /, for their mas|f[ The middle 
panel of Figure |5] shows (for clusters with A/ > 5 x 10^^ A/©) the 
additional fraction of stars that must be formed in this way to ensure 
a consistent model. The discontinuity in the model means that clus- 
ters just above this threshold must form most of their stars in situ, 
since their progenitors all have /»< 0.01. Moreover, most of this 
in situ star formation has to occur at late times, since for most of 
their existence these systems will have had Af < 5 x 10^'^ A/q, and 
they have only recently crossed this threshold. We deal with this in 
the following way. After each merger event, we add sufficient stars, 
via in situ star formation, to bring the cluster back up to the model 
(solid line). No stars are added to (or removed from) systems with 
larger /« than the model. We can then track the average redshift at 
which those additional stars were added, and this is shown in the 
bottom panel of Figure|5] At A/ ~ 5 x Mq, most of the stars 
must be formed very recently, as argued above. This large amount 
of recent star formation in groups is not supported by observations, 
a point on which we will elaborate in ii l3.4| so this is another indi- 
cation that the assumed model is unphysical. For the most massive 
clusters, we only require that ~ 10 per cent of the stars are formed 
in situ, and this at z > 1, which is much more reasonable. 

This fairly complex behaviour nonetheless reflects a consis- 
tent trend in most of our models. The hierarchical merging process 



^ Note therefore the subtle distinction, that a non-evolving /t-A/soo 
model is not synonymous with a lack of star formation. 
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always tends to produce final systems with /* that is nearly inde- 
pendent of halo mass. Thus, for models with a steep h, the most 
massive clusters end up with /, that is too high (and thus incon- 
sistent), while the least massive end up with /, that is much too 
low, and therefore require a lot of recent, in situ star formation. 
This occurs because the mass accretion histories of clusters over 
this limited mass range are not very different - they are built from 
similar-mass haloes over a similar time - resulting in a similar fi- 
nal /, . Thus we can immediately see from this simple example that 
CDM will prefer a value of h that is much closer to than the GZZ 
observations suggest (for non-evolving models). 

We have neglected any discussion of the gaseous component 
of clusters; not only is it dynamically of mino r importance, mak - 
ing up < 20 per cent of the cluster mass (e.g. lAllenet"ai]|2004h . 
but it is expected that maj or mergers are mo st likely to remove gas 
from the dark matter (e.g. lClowe et al.ll2006h . and hence further in- 
crease /*. However, it is also possible th at low-mass clusters and 
groups are relatively deficient in gas (e.g. lAmaud & Evrarjl 19991 : 



Vikhlinin et aL,2006). perhaps because it has been preheated (e.g 



Babul et al.l2002l : lMcCarthv et alj2002l) : if this gas is accreted later 



during the hierarchical growth of structure it could cause /, to de- 
crease. Specifically, if we take the extreme assumption that sys- 
tems with M = 10^^ have no associated gas, then this effect 
alone would lead only to & « —0.04, assuming that systems with 
M — IQ^^ Mq have accreted their full complement of gas. 

This model also implicitly assumes that the radial distribution 
of stars traces that of the dark matter and, in particular, is not altered 
by the merger process. This is a reasonable assumption, as there is 
no evidence that the stellar fraction is a strong function of radius 
outside the very centre of clusters, and the s tellar light distribu - 
tion is usually found to be well modelled by a Navarro et aL Iill997h 



profil e with a reasonable concentration parameter (e.g. 



Lin et al 



2004 ) and comparable to t he total mass distribution jCarlberg et al 



1997 ; iMuzzin et al even for relatively low mass systems 



(e.g. Sheldon et al. 2007). Furthermore, while the merging process 
can greatly distort the relative distribution of gas and dark matter, 
it is much less like ly to separate the stars from the dark matter (e.g. 
IClowe et al.ll2006l) . 



3.4 Constraints 

The previous example demonstrates that there are two aspects of 
our predictions that we can use to choose acceptable models. If 
many systems end up with /, greater than assumed in the model, 
we can confidently rule it out; there is no reasonable way to reduce 
the number of stars, so the model is internally inconsistent. How- 
ever, a model could also be deemed unreasonable if it requires that 
most of the stars are formed too recently. This second constraint is 
less robust because we require guidance from observations. How- 
ever, we will show below that we can afford to be quite conserva- 
tive. 

In Figure [6] we show, for a range of evolution parameters a 
and c, the most negative value of 6 that yields internally consis- 
tent clusters, i.e. those for which the predicted /, on average lies 
on or below the assumed /t-Msoo model. In principle the evolu- 
tion could be positive, such that /, at a given mass scale is greater 
at high redshift than it is today; for example, if efficient star for- 
mation on some scale at 2: = 1 is followed by a long period of 
quiescent accretion of lower mass haloes, in which star formation 
has been inefficient. However, in practice, only negative evolution 
puts interesting constraints on the slope b, so we focus our attention 
on that regime. 




-0.4 -0.6 -0.8 

input, slope b (at. 7. — {]) 

Figure 6. The maximum slope b = d log f*/d log M500 at the present day 
that results in a consistent model with /«< 0.01 in the most massive clus- 
ters. The line colours correspond to the different models shown in Figure[3] 
which make different assumptions about the behaviour of /, for haloes with 
Af < 5 X 10^^ Mq. Dashed lines correspond to evolution in the slope (pa- 
rameter a in Equation]!) while solid lines represent evolution in the normal- 
ization (parameter c). Parameter space to the right of any curve is excluded, 
as it would lead to clusters with too many stars today. The vertical, dotted 
line represents the slope measured from the GZZ data. 



For each model, represented by a curved line, the region to 
the right of the plot is excluded. For non-evolving models, a = 
c — 0, all models show b > —0.33. This means that, given /* = 
0.01 in the most massive clusters, groups with i\f ~ 5 x 10^"^ M0 
must have /« < 0.025. This is considerably shallower than the slope 
found by GZZ, shown as the vertical, dotted line. 

If we allow the 2 = relation to evolve strongly, we can 
weaken our constraints. This is because the groups that merge to 
form clusters will have done so at a higher redshift than the groups 
we are looking at today (recall Figure 0. Figure |6] shows that an 
arbitrarily steep /*-i\f5oo relation at the present day can be accom- 
modated if it is allowed to evolve strongly enough. However, this 
success comes at a high price. In hierarchical models, groups with 
lower mass form at even higher redshifts than the clusters, so they 
will have been built from systems with even lower values of /, . 
The natural prediction of these strongly evolving models is there- 
fore that /, at 2 = will be low not just in clusters, but in all 
haloes. Thus, in order for our model to be consistent, we require 
that a great deal of stars form in situ in these lower mass haloes, af- 
ter most of the clusters have been assembled, which, for a ACDM 
Universe, is 2 ~ 0.3 (see Figure |9}. The more steeply we assume 
/, evolves, the more recently those stars must have been created. 

We quantify this in Figure [7] For each halo at 2 = 0, we as- 
sume that any stars that were generated via in situ star formation 
formed at the first redshift that any merger yields a stellar fraction 
below the assumed valu^. Stars that were accreted through merg- 
ers are conservatively assumed to have formed at 2 — 2. We use 



* In reality, the stars could form anytime between then and the next merger 
This means our results are conservative. 
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Figure 7. The median redshift at which stars must have formed in haloes 
with 5 X 10^"^ < M/Mq < 10^*, to remain consistent with each model, 
as a function of the maximum allowable slope b at the present day, from 
Figure |6] The line styles correspond to those in Figure |6] For each model, 
more negative values of b imply greater evolution in the /»-Af5oo relation. 
The vertical, dotted line represents the slope measured from the GZZ data. 



this to compute, for each halo, the redshift at which half the stars 
were formed in situ, zsf. Figure|7]shows the median zsf for haloes 
with 5 X 10^^ < M/Mq < 10^" (i.e. at the low mass end of our 
models), as a function of the maximum z = slope b allowed 
by each model, from Figure (6] Recall that more negative values 
of b imply steeper evolution. As described above, the success of 
strongly-evolving models to match the low /, in today's clusters 
comes at the expense of requiring substantial recent star formation 
in groups. 

Although it is clear that galaxy cluster pop ulations are 
old and passively evolving, with zf S> 2 (e.g. iBower et al 



1993; 



2001 



De Propris et"al]|l999l, [20071; iFinn et al] |2005| ; iNelan et al 



Lin etalj|2003 ;l Muzzin et al. I l2007l ). there are fewer con 



straints on lower mass groups. However, globally we know that 
the total mass in stars has at most doubled since z = 1 



(e.e.. iDickinson et all |2003|; iBell et alj 


20031; iDrorv et all |2005|; 


iGwvn & Hartwickl 20051; Pozzetti et al. 


2007h. Relative to the 



global population, at z < 0.5 groups are known to be dominated 
by galaxies w ith early-morph ological typ es and little or no star for- 
mation (e.g. jde la Rosa et al. 2001 ; Tran et al. 2001 ; Balogh et al. 
20021 l2004l2007.; .Weinmann et al. ,2006. ; .Wilman et aL ,2005 ; 



Jeltema et al l2007l ; lMcGee et al.l2007h ; herefore we would expect 
them to h ave formed half their stars well before z — 1. Recently, 
iBrough et al. (2007) have made a detailed analysis of the brightest 
galaxies (which, together with the ICL, GZZ claim dominate the 
stellar mass) in three X-ray groups. Two of these have luminosity- 
weighted ages > 10.5 Gyr (la) limit, corresponding to a formation 
redshift of z/ = 2; the youngest has an age limit of > 6.9 Gyr, or 
Zf = 0.75. 

We therefore consider that zsf ^ 1 (a lookback time of 8 Gyr) 
is a very reasonable lower limit for the redshift at which groups 
(5 X 10^^ < M/Mq < lO") have formed most of their stars. 
From Figure[7]this implies a lower limit of fe > —0.35, only slightly 



steeper than the non-evolving model constraints. To accommodate 
a slope b — —0.64, as observed by GZZ, would require that at 
least half the stars in groups with 5 x lO" < M/Mq < 10^'* 



formed since z 
observations. 



0.35, which is very unlikely given the above 



3.5 Ab-initio models 

Instead of assuming a correlation between /, and M500, an alter- 
native approach is to investigate more complex, ab-initio models 
that include prescriptions for star formation and feedback. There 
are many su ch models currently available, with generally similar 
recipe s (e.g. lOe Lucia et al.ll20o3 ; ICroton et alfcOOdlBower et all 
20061). Here we consider the publicly availabl4j predictions of the 



Bower et al. I ( l2006h model. Figure [8] shows the predicted /* dis- 
tribution as a function of halo mass, for all systems in the par- 
ent simulation, at redshifts between z = and z = 2. For halo 
mass, we use a friends-of-friends mass with linking length b — 0.2, 
which corresponds approximately to -M200, with a scatter of about 
15 per cent. For a typical cluster halo, M5oo~ 0.65M2oo- There 
is almost no mass dependence of /, for A/ > 10^'^ M0, and very 
little evolution in the relation. Above M — 10^'^ M0, the model at 
all redshifts satisfies approximately d log log Afsoo ~ —0.05. 
This prediction for a nearly constant /, is also in good agreement 
with cosmological hydrodynamic simulations that includ e coo ling 
and feedback physics (e.g. Borg ani et alT2004) ; lKav etliLl2007t) . al- 
though those simulations tend to predict a higher overall value of 

This prediction for a nearly constant /, on these scales is not 
surprising, given that massive groups and clusters are both built 
from haloes with a similar mass distribution, and over a simi- 
lar timescale. Although our arguments from the previous section 
demonstrate that CDM could support a ./*-M5oo relation as steep 
as 6 ~ —0.35 at the present day, this is only true under the most 
conservative, and probably unrealistic, assumptions about the value 
of /* in haloes with Af5oo< 5 x lO^'^AfQ and how this evolves with 
redshift. 



3.6 Other cosmologies 

We can hope to put constraints on the cold dark matter model in 
general, regardless of the specific cosmology, since they are all hier- 
archical in nature. Different cosmological parameters will result in 
different rates of formation for haloes of a given mass. We showed 
in i; l3.4l that one can maintain a low /» in the most massive clusters 
today, whatever the local /.-A/soo relation, as long as this relation 
evolves strongly with redshift. Clearly one could achieve this same 
result with weaker evolution, in a cosmology where structure forms 
earli er. This is true of an o pen Universe, with Qrn < 1 and f^A ~ 
(e.g. Ivan den Bosclill2002h ; alternatively the epoch of structure for- 
mation can be pushed to higher redshif t if th e normalization of the 
power spectrum, as, is increased (e.g. iLacev & Cole 199!%_This 
is illustrated in Figure |9] where we use the iLacev & Cold l ll993h 
model to calculate the probability distribution of the formation red- 
shift, defined as the redshift where 75 per cent of the final mass is 
in place (see details in Balogh et al. 1999). Distributions are shown 
for haloes with final masses of lO^Afg, IO^^Mq, and 10^^ Mq, 
for a range of cosmological parameters. As expected, the highest 
formation redshifts are obtained in a low-Qrn or high-as Universe. 
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Figure 8. The predicted stellar fraction as a function o f halo mass (approx - 
imately M200), from the galaxy formation models of iBower et alj j2006l) . 
The points represent binned averages at different redshifts: blue, green, yel- 
low, and red points correspond to z = 0, 0.5, 1.0, and 2.0, respectively. 
The blue shaded region corresponds to the 10th and 90th percentiles of the 
z = distribution. The dashed line is an approximate representation of the 
z < 1.0 model data for Afsoo > lO^^M©, and has a slope of b = —0.05. 
The dotted line represents the approximate slope (b = —0.64) of the GZZ 
data, with arbitrary normalization. 



However, in all cases, lower-mass haloes on average form at 
even higher redshifts; this of course is the well-known behaviour of 
cold dark matter models. Therefore, in a low-f^m or high-as Uni- 
verse, galaxy groups had most of their mass in place at even higher 
redshift than the clusters, and /« today should also be lower, re- 
quiring considerable late-epoch star formation in groups to retain 
consistency with a steep relation b <^ 0. This trade-off means that 
our constraint on b is unlikely to be changed in different cosmolo- 
gies. 

We have rerun our Pinocchio simulations, for an Einstein De- 
Sitter Universe {Q.m ~ 1, S7a = 0, ag = 0.5) and an open Uni- 
verse (flm = 0.1, flA = 0, (T8 — 0.9). To be conservative we 
adopt the Min model for the f^-M^oo relation at 2; = 0, and allow 
the normalization or slope to evolve rapidly, as (1 + z)^. The min- 
imum value of b in this model, that is consistent with /,= 0.01 in 
the most massive systems and ensures that at least half of the stars 
formed at 2 > 1, is still b ~ —0.3 for the EdS model, but could 
be relaxed to 6 ~ —0.4 for the OCDM cosmology, only slightly 
steeper than the constraint for our default ACDM case. Thus, our 
conclusions are nearly independent of cosmological parameters, as 
expected. 



4 DISCUSSION 

4.1 Observational Uncertainties 

As published, the GZZ data show stellar fractions reaching as high 
as 30% in the lowest mass systems (which, we note, have poorly 
calibrated M500 ), well above the WMAP c onstraints on the global 
baryon fraction dSpergel et al]|2003ll2007h . If this is confirmed to 
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Figure 9. The probability distribution of formation redshifts, defined as the 
redshift where 75 per cent of the mass is assembled, for clusters with virial 
masses M = IO^^Mq (sohd line), M = (dotted line), and 

M = 1 0^^ Mp, (dashed line) . The calculations are based on the analytic 
model of'Lacev & Coi3 il993l) . Different panels show the results for differ- 
ent cosmological parameters, as indicated. The highest formation redshifts 
are achieved in low-Hm or high-erg models. But in all cases, low-mass 
haloes form at substantially higher redshift than high-mass haloes. Thus, 
cosmological parameters have little effect on the generic prediction of hier- 
archical models, that b = d log /* /rf log A/500 as measured at the present 
day must be flat, b > —0.35. 



be representative of systems in this mass range, it will rule out hi- 
erarchical structure formation models. The shallower relation be- 
tween /* and A/500 as found by LM, however, does appear to be 
consistent with such models. In this section we will look carefully 
at some of the biases and uncertainties involved in the analysis. 

One of the most important differences between these two sur- 
veys is in the measurement of total mass. GZZ derive masses from 
the velocity dispersions, which have significant uncertainties them- 
selves, and translate into relative errors three times larger when 
converted to M^oo- Moreover, the strong correlation between er- 
rors on A/500 and Ai^,,5oo, as discussed in §|2j2] further complicates 
matters. To illustrate this better we replot the data in Figure[TO] this 
time showing A/500 as a function of A/*, 500, so that the tilt in the 
error bar now reflects the degree of correlation between these two 
measurements. 

It is evident that the error bars have the same "slope" as the 
GZZ data itself, suggesting that the steepness of their relation is 
at least partly due to this correlation. However, this may not be 
the whole story. Accounting for the correlation on the errors, we 
can compute how likely it would be to find as many systems with 
/»> 0.05 as GZZ do (seven), if they all actually had /.= 0.02 
(consistent with LM), and random uncertainties on a scatter the 
observations. This probability is less than 0.1 per cent. The uncer- 
tainties on A'/5oo do not include the scatter in the a-M^oo relation 
from which it is derived; however, even increasing the error bars 
on a by 50 per cent, we find that we would only expect to find 
seven systems with /*> 0.05 six per cent of the time. It is clear 
from Figure[TO]why this is; the GZZ data lie systematically below 
those of LM; while any one point may be discrepant by only one or 
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Figure 10. The total mass within 500 is shown as a function of the total 
stellar mass (including intracluster light) within the same radius, for the LM 
sample (squares) and the GZZ sample (circles). The two open circles rep- 
resent clusters A2405 and APMC020. which are systems strongly affected 
by line-of-sight structure. The dashed lines show lines of constant stellar 
fraction, as indicated. Error bars reflect the uncertainty on M500, and the 
tilt shows how this is correlated with A/».5oo. 



two standard deviations, the difference between the two relations is 
much more significant. If the source of the discrepancy lies in the 
total masses, then, it implies that these masses are systematically 
underestimated for Msoo;^ 10^''. One possibility is that the veloc- 
ity dispersions themselves are underestimated; however, this seems 
unlikely, as several authors have foun d that <j is quite robust for sys- 
tems with at least 20 members (e.g. Zabludoff & Mulchaevl[T99g : 
iBorgani et alll 19991 : iHicks et alj|2006h . The other, more likely ex- 
planation is that the adopted calibration between cr and Mr,oo is 
incorrect at low masses. In particular, the four groups with the low- 
est a require an extrapolation of this relationship, which may well 
be unreliable. Ignoring these four clusters, the remaining GZZ data 
are statistically consistent with a stellar fraction of 0.02 that is in- 
dependent of cluster mass, and an apparently steep slope that is due 
entirely to the correlated error bars. 

There are other possible sources of systematic error in the 
data, but none of them seem likely to account for the discrepancy. 
The most obvious place to look is the intracluster light, since this 
was not directly measured by LM. In fact, the claim by GZZ that the 
ICL fraction is such a strong functio n of mass is not rea dily appar- 
ent in the recent measurements from lZibetti et all ( l2005h . based on 
ensemble averages of SDSS clusters. Moreover, numerical simula- 
tions generally find that the ICL component s hould actually be less 
impo rtant in groups, relative to clusters (e.g. iMurante et alj|20()4l 
l2007h . However, even though GZZ claim the ICL doubles the stel- 
lar mass in the smallest observed systems, this is still not enough 
to account for the discrepancy with LM. In Figure[TT]we replot the 
data shown in Figure [TO] but excluding the ICL component from 
both samples. There is still a significant discrepancy between them, 
and the lowest mass systems in GZZ have stellar fractions > 10%. 
In particular, if the ICL is ignored, then LM observe a nearly con- 
stant /, over all masses, not only consistent with our theoretical 



M. 



(excluding ICL) 



Figure 11. As Figure [Tol but where we have excluded the ICL contribution 
to the stellar light, for both surveys. 



bounds but also in good qu antitative agreement with ab-initio mod- 
els (e.g. lBower et alfcood) . On the other hand, GZZ still predict a 
strong mass dependence of the fractional galaxy light alone. 

The dominant systematic uncertainty in the stellar mass mea- 
surements are likely in the stellar mass-to-light ratios. The abso- 
lute, average, M/Li adopted is of little consequence, but the rela- 
tive difference between M/Li and M/ Lk, and any trend in these 
values with total mass, is relevant. The relative value of M/Li 
and M/Lk, and hence the relative normalization of the data in 
Figs |2l [To] and [TT] depends on our assumption that the average 
galaxy colour is I ~ K — 2.0 in these samples. A bluer pop- 
ulation would act to decrease the GZZ stellar masses, relative to 
LM, and thus decrease the systematic offset observed in these fig- 
ures. However, no reasonable colour will reconcile the four lowest- 
mass GZZ groups with those of LM, and the remainder are already 
consistent with LM (and a constant stellar fraction) given the er- 
ror bars. Of more interest is the possibility that M/L varies with 
system mass. If lower-mass clusters have systematically younger 
stellar populations they will have lower average M/L, and if un- 
accounted for this would lead to an artificially steep slope b (and 
would also compromise some of the conclusions in GZZ). How- 
ever, to explain the discrepancy with the LM data would require that 
the low-mass clusters of GZZ are systematically bluer than those of 
LM, while the massive clusters are similar. It seems unlikely that 
either of these biasses are important, since GZZ find that most of 
the stars in these low-mass systems are in the ICL and BCG com- 
ponent, and the BCG at least i s not likely to be mu ch younger than 
in more massive clusters (e.g. iBrough et ai]|2007h . Moreover, the 
average stellar M/L in the / band is unlikely to vary by more than 
a factor of about two, which is insufficient to reconcile the steep 
slope b observed by GZZ with the model predictions. Nonetheless, 
deep infrared images of these systems would be very useful. 

Another concern might be the contribution to the stellar light 
from faint, unresolved galaxies; however, both GZZ and LM as- 
sume fairly steep, mass-independent faint end slopes when extrap- 
olating the galaxy luminosity function (a = —1.21 and a = —1.1, 
respectively), so this cannot contribute to the difference. Finally, 
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there is a possible selection effect as GZZ select groups to have 
a dominant galaxy, and this may bias them toward systems with 
particularly high stellar fractions. Observations of a more repre- 
sentative sample would be valuable, especially if there are suffi- 
cie nt redshifts to rob ustly identify cluster members. However, in 
the lBower et alj ( 1200 6) models we considered in §[33] none of the 
systems in the relevant mass range are found to have stellar frac- 
tions as high as even 5%. If these models provide an accurate pic- 
ture of the local Universe in this respect, then no possible selection 
bias would lead to the high values of /, observed by GZZ. 

4.2 The consequences 

If the strong trend observed by GZZ is confirmed, and found to 
be typical of a mass-limited sample of groups, what are the theo- 
retical consequences? We have shown that such groups cannot be 
the progenitors of today's clusters, even if they have built up their 
high stellar fraction quite recently. One possible implication would 
be that today's clusters have grown very little in mass since at least 
z = 1, while groups have assembled more recently. This inherently 
non-hierarchical model would require a lot of suppression of power 
on group scales; qualitatively this is the behaviour of warm dark 
matter models, alth ough these generally suppress structure on much 
smaller scales (e.g. lAvila-Reese et alJl200ll : lBode et alj200ih . An- 
other possibility is that a large fraction of massive dark matter 
haloes have no stars associated with them, or at least not enough 
to be det ectable. There is some preliminary evidence for su ch dark 
clusters dvon der Linden et al. I l2006l ;l Mahdavi et al. ||2003). How- 
ever, for this to be the solution, the most massive clusters would 
have to have accreted ~ 80% of their mass from such objects, while 
groups still accrete all their mass in haloes with /, > 0. This seems 
quite unlikely, but should be testable in forthcoming weak-lensing 
surveys. 

On the other hand, if the dynamical masses of the groups in 
the GZZ sample are underestimated, this raises another theoreti- 
cal problem. A particularly nice result from GZZ is that (ignoring 
the four lowest mass clusters) their data give a full account of the 
baryons: the sum of the stellar mass and the expected gas mass 
(unfortunately, not directly measured) is comparable to the total 
baryon fraction in the Universe, and independent of system mass. 
This is an attractive explanation for the observation that galaxy 
groups are deficient in X-ray gas (e.g. lAmaud & Evrard Il999l : 
[vikhlinin et al. 2006): the missing gas has cooled to form stars. 
If GZZ have overestimated however, this explanation will no 
longer be viable, and the low gas fractions of groups would likely 
imply that some powerful form of heating has expelled a large frac- 
tion of the gas beyond 500. However, this heating must not be too 
strong, so that the most massive systems are still able to retain all 
their gas. This requires something of a delicate balance. 



5 CONCLUSIONS 

The stellar fraction /, as a function of cluster mass Afsoo is 
an important test of hierarchical structure formation models. We 
find that such models can make a robust, falsifiable predic- 
tion that the power-law slope relating these two quantities is 
d log /*/ci log i\f5oo > —0.35, for systems with A/soo > 5 x 
1O^''M0. a steeper slope can only be accommodated if the 
M500 relation evolves strongly, such that galaxy groups formed 
most of their stars in situ since z = 1.0, which is not sup- 
ported by observations. Since the most massive clusters today (with 



~ 10^^ Mq) are robustly measured to have /,~ 0.01, 
hierarchical models therefore require that galaxy groups (with 
M500 ~ 5 X 10^^ Mq) have /*< 0.03. This constraint is a con- 
servative limit; ab-initio models pr edict a much flatter relationship 

with h > -0.1 teower et alj2006h. 

Recent observations bvlO onzalez et alj j2007h appear to con- 
flict with this prediction. In particular, their data follow a power- 
law relation over almost two orders of magnitude in mass, with 
6 — —0.64, and the lowest-mass systems in their sample have stel- 
lar fractions of 30 per cent, exceeding even the limits on the baryon 
fraction from WMAR If confirmed, these observations definitively 
rule out hierarchi c al stru cture formation. A"— band observations 
from iLin &Mohil j2004h . on the other hand, are just consistent 
with the model co nstraints, but may s till be inconsistent with the 
ab-initio model of iBower et alj ( |2006|) if the ICL contrib ution is 
as dominant in groups as claimed bvl Gonzalez et alj j2007h . These 
data do not extend to such low mass systems, nor are their data deep 
enough to directly measure the important ICL contribution. More 
observations are needed to resolve the discrepancy between these 
two studies, and to thereby test the viability of cold dark matter 
models. Particularly valuable would be X-ray images (and temper- 
ature maps) of the Gonzalez et al. ( 2007) sample, to measure the 
gas content and total m ass, and deep / or near-infrared images of 
the lLin&MohiH2004) systems to directly measure the ICL com- 
ponent. 

We end by noting however that the four lowest-mass systems 
in the Gonzalez et al. sample have poorly calibrated masses. If we 
accept that their masses are underestimated by a factor of ten, then 
the remainder of their data are consistent with a universal stellar 
fraction of a few percent, and the apparent correlation with system 
mass is due to the correlated error bars on dynamical and stellar 
masses. Furthermore, removing these four systems greatly reduces 
the significance of their claim that the intracl uster light fraction is 
a strong function of mass and, in this case, the lLin & MotiJ j2004h 
data (which have much smaller error bars) are also consistent with 
a constant stellar fraction /,~ 0.01. This conclusion challenges 
the claim by GZZ that the stellar component of galaxy groups is 
dominated by the BCG and ICL, and that the low gas fractions in 
these groups is attributable to an increased stellar fraction. 
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